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Abstract 

We show that high energy scattering is a statistical process essentially similar to 
reaction-diffusion in a system made of a finite number of particles. The Balitsky- 
JIMWLK equations correspond to the time evolution law for the particle density. 
The squared strong coupling constant plays the role of the minimum particle den- 
sity. Discreteness is related to the finite number of partons one may observe in a 
given event and has a sizeable effect on physical observables. Using general tools 
developed recently in statistical physics, we derive the universal terms in the rapid- 
ity dependence of the saturation scale and the scaling form of the amplitude, which 
come as the leading terms in a large rapidity and small coupling expansion. 



1 Introduction 



Much progress has been made recently in understanding high energy hard scat- 
tering in QCD at or near the unitarity limit. General equations have been given 
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by Balitsky [1] and by Jalilian-Marian, Iancu, McLerran, Leonidov, Kovner 
and Weigert (JIMWLK) [2,3,4] which generalize the Balitsky- Fadin-Kuraev- 
Lipatov (BFKL) evolution [5] to the region where unitarity (saturation) ef- 
fects become important. The Balitsky- JIMWLK equations are nonlinear op- 
erator equations, while a "mean field" version of the equations, the Balitsky- 
Kovchegov (BK) [1,6] equation, is a nonlinear equation for the scattering am- 
plitude and has been widely studied recently. The scattering amplitude which 
emerges from the BK equation is, in general terms, characterized by the energy 
(rapidity) dependence of the saturation momentum Q S (Y) [7,8,9,10,11,12], and 
by geometric scaling [13,9,10], the statement that the scattering amplitude 
A(Q 2 , Y) is equal to a function of a single variable A(Q 2 /Q 2 (Y)). 

However, it has not been clear to what level the general properties of solutions 
to the BK equation are shared by solutions to the Balitsky- JIMWLK equa- 
tions. This is the problem we address in this note. Our object is to describe 
the energy-dependence of the saturation momentum and the scaling properties 
of the scattering amplitude which should emerge from the Balitsky- JIMWLK 
equations. When viewed in a particular way the problem here looks identi- 
cal to a class of problems studied recently in statistical physics [14,15]. In 
the statistical physics problems the Fisher-Kolmogorov-Petrovsky-Piscounov 
(FKPP) equation [16] approximately describes the time evolution of certain 
quantities in some discrete systems. In the QCD problem of dipole-dipole 
scattering the BK equation describes the rapidity and Q 2 dependence of the 
scattering amplitude. As has been noticed recently, the BK equation is in the 
same universality class as the FKPP equation and this fact gives a powerful 
and general derivation of the energy dependence of Q S (Y) and of geometric 
scaling [12]. 

The FKPP equation has limitations in applications to average quantities in 
discrete statistical systems. When the discreteness is not important, that is in 
a region where many "objects" are present, the FKPP equation is a good ap- 
proximation to the actual evolution of the system. However, when only a few 
objects are involved discreteness effects are significant and the FKPP descrip- 
tion breaks down [14,17]. A similar effect occurs in QCD evolution, and this is 
seen most easily by viewing the scattering of an elementary dipole of size r on 
an evolved dipole of initial size r in terms of the Y — evolution of particular 
configurations of dipoles starting from r®. (The scattering amplitude is then 
given by an average over all possible configurations.) The dipoles making up a 
configuration are the discrete elements of our system. The importance of fluc- 
tuations due to discretness in QCD evolution was first noticed by Salam from 
Monte-Carlo studies [18]. More recently, the importance of fluctuations has 
been reiterated in the context of non-linear evolution in Ref. [19], where the 
role of rare fluctuations in the approach of the S-matrix towards the unitarity 
limit has been discussed, and also in Ref. [4,20], where JIMWLK evolution 
has been reformulated as a random walk in some functional space, thus em- 
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phasizing its stochastic nature. This last formulation lies also at the basis of 
the numerical studies of JIMWLK evolution in Ref. [21]. 

So long as the dipole occupancy in a configuration is large compared to one, 
use of the BK equation should be a good approximation for the evolution of 
our configuration. However, when there are only a few dipoles of size r in 
our particular configuration BK evolution cannot be expected to be accurate. 
Indeed in a discrete picture occupancy can go below one only by becoming 
zero which stops the evolution along that path. Thus, using the BK equation 
with a cutoff when dipole occupancy is near one, or equivalently when the 
scattering amplitude for a particular configuration becomes of size a 2 s , exactly 
the same procedure used for discrete statistical systems [14,17], should be a 
good representation of the evolution of the system. This cutoff is essentially 
the same as that introduced in Ref. [22], and the present discussion can be 
viewed as a justification of the procedure used there at least for the calculation 
of the energy dependence of Q S (Y). 

To also compute the dependence of the scattering amplitude upon the dipole 
size, and thus compare with the geometrical scaling form of the solution to 
BK equation, one needs to understand the fluctuations of the saturation mo- 
mentum from one configuration to another. Here our control is less complete, 
and we rely on a scaling law recently seen in numerical simulations. The scale 
which emerges is equal to the square root of the value found in Ref. [22] where 
fluctuations at the boundary were not included. 

Finally, it should be emphasized that our description is for a scattering at 
a definite impact parameter. A more complete discussion, including impact 
parameter dependences, will be given later [23]. 



2 High energy scattering as a statistical process 

We consider the scattering of a dipole of variable size r (the probe) off a 
dipole of size r (the target). A natural variable that will be used throughout 
is p = ln(rg/r 2 ). We go to the rest frame of the probe so that the target carries 
all the available rapidity Y . The impact parameter b between the dipoles is 
fixed. 

The target interacts through its quantum fluctuations, which at high energy 
are dominated by gluons. It proves useful to represent this set of partons by 
color dipoles [24]. This is possible in the dilute regime in which saturation 
effects (i.e. interactions among gluons inside the target wavefunction) are neg- 
ligible, but the effects of fluctuations should on the contrary be important. 
Then, the dipole picture emerges in the large- N c limit, in which gluons are 
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similar to zero-size qq pairs and non-planar diagrams are suppressed. The 
dipole approximation breaks down when the amplitudes approach their uni- 
tarity limits (indeed, in the considered frame, the unitarity corrections are 
tantamount to saturation effects in the target). However, that does not ham- 
per getting the right asymptotics for physical quantities like the saturation 
scale since, as we shall see, this is controlled by the dynamics in the tail of the 
distribution at high transverse momenta, or small dipole sizes. 

We denote by T(r, r ) the scattering amplitude of the probe off a given par- 
tonic realization \uj) of the target (the dependence on b is understood). It is a 
random variable, whose probability distribution is related to the distribution 
of the different Fock state realizations of the target. The values of T(r) range 
between (weak interaction) and 1 (unitarity limit). T(r) will be an essential 
intermediate quantity in our calculations, but it is not an observable. The phys- 
ical dipole-dipole scattering amplitude A(r, Y) is the statistical average over 
all partonic fluctuations of the target at rapidity Y, i.e. A(r,Y) = (T(r))y. 

When T is small, T(r, r ) = T e i{r, rj), where % labels the dipoles in the Fock 
state of the target at the time of the interaction. T e i is the elementary dipole 
interaction and is essentially local in impact parameter. T d behaves like 

T e i(r,ri) ~ a 2 s , with r< = min(|r|, |r;|) , r> = max(|r|, |r;|) (I) 

when the dipoles overlap, and vanishes otherwise (see the insert in Fig. I). 
We have neglected 0(1) factors and logarithms, but these approximations do 
not affect the results that we shall obtain, which are largely independent of 
the details. Eq. (1) shows that the amplitude T(r, r ) is simply counting the 
number n(r, r ) of dipoles of size r within a disk of radius r centered at the 
impact parameter of the external dipole (the dipole occupation number): 

T(r,r ) ~ a 2 s n(r,r ) . (2) 

Note that n(r, r ) can take only discrete values. Thus, in this description, 
fluctuations in T emerge naturally as fluctuations in the particle (here, dipole) 
number, which should be especially important in the regime where n ~ (9(1). 
The unitarity bound on T implies that n(r, ro) is also constrained by an upper 
bound N ~ 0(1/ a%). 

The dipole picture also provides us with the evolution law for the dipole dis- 
tribution with increasing rapidity. Consider a small increment dY of the total 
rapidity from a boost of the target. Then, each of the dipoles r already present 
in the wave function from the previous evolution and for which n(r, r ) < 1/ a 2 s 
may split into two new dipoles, of respectives sizes z and r — z, with a differ- 
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ential probability [24] 



a 



r 



2 



dP = dY 



2n z 2 {r- — z) 



-^d 2 z = [XadY] x \p(z, r — z\r) d 2 z] 



(3) 



where a = a s N c /n. We see that aY is the natural evolution variable: we will 
call it "time" . The second equality in Eq. (3) expresses dP as the product of 
the inclusive probability of splitting 3 



in the time interval udY, by the conditional distribution of the sizes of the 
produced dipoles 



By applying Eqs. (3)-(5) to any of the dipoles r\ present in the wavefunc- 
tion, one can easily deduce the evolution law for the dipole configuration as a 
whole. This leads to a description of the (dilute tail of the) target wavefunction 
as a stochastic ensemble of dipole configurations endowed with a probability 
distribution which evolves with Y according to a master equation [25]. This 
picture, which is similar to certain problems in statistical physics, is partic- 
ularly appropriate for a study of fluctuations, since the discreteness of the 
dipole number is explicit. However, this picture breaks down, as anticipated, 
when the dipole occupation numbers — which in the dilute regime rise expo- 
nentially with Y — become of 0(l/a 2 ), and saturation effects start to play a 
role. But in this high density regime, one can rely on a different formalism, the 
color glass condensate [3], which is an effective theory for gluon correlations in 
the target wavefunction at small x, and is endowed with a functional evolution 
equation — the JIMWLK equation [2,3,4] — which shows how these correla- 
tions change under a boost. When this formalism is applied to the scattering 
between the color glass (the target) and a set of dipoles (the projectile), the 
saturation effects encoded in the JIMWLK equation are converted into uni- 
tarity effects in the evolution of the scattering amplitudes. The latter are thus 
found to obey an infinite hierarchy of evolution equations originally derived 
by Balitsky [1]. 

In what follows, we shall need only the first equation in this hierarchy, which 
applies when the projectile is a single dipole. This equation is most easily ob- 
tained by using the rapidity increment dY to accelerate the projectile (the 
dipole of size r). Within the rapidity interval dY, either the dipole does 
not split, in which case its scattering amplitude T(r) = T(r, r ) remains un- 
changed, or it splits, in which case T(r) is replaced by the scattering amplitude 

3 An ultraviolet cutoff is understood in Eq. (4). It disappears in physical quantities. 
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p(z, r — z\r) d 2 z = A 
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Fig. 1. The amplitude T for a typical partonic realization as a function of 
p = ln(rg/r 2 ). The individual dipoles seen at impact parameter b are represented 
by a short vertical line. The straight line is the sum of their contributions to the 
amplitude. In the saturation regime, the dipole description breaks down, that is 
indicated by the filled box. Upper right corner: the contribution of a single dipole 
to T, Eq. (1). 

of the two child dipoles. This leads to the following evolution law 

j T{t)\y with probability 1 — A adY 

{r)\ Y +dY = +T ( r _ z ) _ T ( z ) T ( r _ z )| y with probability XadY 

(6) 

where z is distributed according to p(z,r — z\r) d 2 z. Taking the limit dY — > 
and replacing A and p from Eqs.(4),(5), one gets 4 

dr(T(r)) Y = £ / d*z^-y 2 ((T(z)) Y + (T(r-z)) Y - (T(r)) Y 

-(T(z)T(r-z)) Y ). (7) 



As anticipated, Eq. (7) is not a closed equation for (T): it depends upon 
the correlator (T(z)T(r — z))Y- A mean field approximation (T(z)T(r — z)) ~ 
(T(z))(T(r— z)) would cast Eq. (7) into a closed form, known as the Balitsky- 
Kovchegov (BK) equation [1,6]. The linearized form of Eq. (7) is recognized 
as the (dipole version of) BFKL equation [5]. 

Let us finally discuss the typical shape of T{r) as resulting from the previous 
considerations. It is a well known characteristic of the BFKL evolution that 



4 The impact parameter dependence could be easily put back in Eq. (7). We have 
omitted it for simplicity and since it is enough for our purpose to assume locality 
of the evolution. 
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the most probable splittings are those in which both child dipoles have a 
size comparable to the size of the parent dipole (this can be also checked on 
Eq. (7)). Thus, if one starts with one dipole r at Y = 0, then the main 
mechanism for the rise of T(r) with Y is a growing diffusion around the size 
of the initial dipole r . Consequently, in a typical partonic configuration as 
obtained after a sufficiently large rapidity evolution, the dipoles appear to be 
densely distributed around the size r$ (where T(r) is large), but they become 
more rare with decreasing r (or increasing p), and for sufficiently large p one 
meets only rare fluctuations which involve one (or few) dipoles and for which 
T(r) ~ a 2 . The typical partonic realization is shown in Fig. 1: it is a front 
which with increasing Y progresses towards larger values of p. We define the 
saturation scale Q S (Y) of a given partonic configuration by the position of 
this front, that is, by the value of the inverse dipole size for which T reaches 
some predefined number T of order one: T(1/Q S (Y)) = T . We also define 
Ps(Y) = \og{rlQl{Y)). 



3 The energy dependence of the saturation scale 



When T(r) is of order a 2 , the number of dipoles participating in the scattering 
is small, and fluctuations dominate the dynamics of T(r). By contrast, when 
T(r) ^> a 2 , the fluctuations 5T in T become relatively unimportant (since 
typically 5T ~ a s VT), so the dynamics is self-averaging (for each individual 
front realization), and the mean field description of a given event becomes 
justified. Consequently, the evolution of T in the bulk of the front T(r) 3> a 2 
is essentially given by a mean field equation (the BK equation). It turns out 
that, for the purpose of computing the asymptotic energy dependence of the 
saturation momentum, one can still rely on a modified mean field approxima- 
tion, which is obtained by introducing a factor 0(T — a 2 ) in the BK equation. 
The latter implements the fact that, in a real event, in which the occupation 
number is discrete, T cannot become less than a 2 (cf. the discussion in Sec. 2). 

To appreciate the dynamical role of this cutoff, it is useful to notice an essential 
difference between the tail of a real event, and that of the solution to the BK 
equation: whereas the front generated by the BK equation has an exponential 
tail which extends up to arbitrarily large p (see Eq. (8) below), a real event, 
on the other hand, is like a histogram whose front is necessarily compact: for 
any Y, there exists a foremost occupied bin (f.o.b.) pf. G .b. = Pf.o.b.(^) suc h 
that T(pf.o.b.)|y > and T(p)\ Y = for any p > pf. Q .b.- This implies that the 
mechanism for front propagation is different in the two cases. For the mean 
field approximation, the dominant mechanism is the local growth within the 
tail of the distribution: at any p ^> p s (Y), the local amplitude rises very fast 
due to the BFKL instability, thus "pulling" the front towards the right. By 
contrast, in a real event, the local growth is not possible in the empty bins on 
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the right of the f.o.b., so the only way for the front to progress into those bins 
is via diffusion, i.e. via radiation from the occupied bins at p < Pf.o.b.- But 
since diffusion is less effective than the local growth, we expect the "velocity" 
dp s (Y)/dY of the front (i.e., the exponential growth rate for the saturation 
momentum) to be reduced in the real event as compared to the corresponding 
prediction of the BK equation. 

A simple way to try and capture this physical situation in mathematical terms 
is to insert a cutoff 0(T — a 2 s ) on the growth term in the BK equation, while 
allowing the diffusion there to remain operative even at arbitrarily small T. 
(This is possible after separating the local growth term from the diffusion term 
in the BFKL kernel with the help of a "diffusion approximation"; see, e.g., [12] 
for details.) This simple recipe was in fact invented by Brunet and Derrida [14] 
in the context of statistical physics. They studied the propagation of fronts in 
the presence of fluctuations associated with discreteness in a variety of physical 
situations (see Ref. [17] for a review). Using the mean field approximation to 
their dynamical equations together with a suitable cutoff, they were able to 
compute analytically the time evolution of the position of the front and of its 
bulk shape. Although there is so far no full mathematical justification of this 
procedure, it has been checked (through numerical calculations) that it yields 
indeed the right value for the velocity of the front at large times and for a 
large number of particles. This result is likely to be valid for all models that 
fall in the universality class of the stochastic FKPP equation [17]. 

We expect that, also in the present QCD context, that recipe give the right 
asymptotics of p s (Y) = ln(rlQ^(Y)) for ln(l/c^) 1. Indeed, when viewed in 
the way exposed in Sec. 2, the rapidity evolution of T is essentially the same 
as the time evolution of say the particle number density of a system made of 
a number N ~ l/c^ of diffusing and interacting particles, and hence belongs 
to the class of models studied by Brunet and Derrida. 

We now turn to the practical computation of the rapidity dependence of the 
saturation scale. Since it amounts to solving the BK equation supplemented 
by a cutoff, it is useful to recall first the asymptotic solutions to the BK 
equation without a cutoff and the way how they set in (see e.g. [12] for 
details). As the BK equation falls into the universality class of the Fisher- 
Kolmogorov-Petrovsky-Piscounov (FKPP) equation [16], its asymptotic solu- 
tion is a traveling wave, i.e. a uniformly translating front whose "position" is 
characterized by the logarithm of the saturation scale p s (Y). It moves with 
the "velocity" dp s (Y)/dY = ax(lo)/lo towards smaller dipole sizes. Here, 
x{l) = 2-^(1) — VKt) ~~ "^(l - 7) 1S the p-moment of the dipole splitting prob- 
ability (3), or, equivalently, the characteristic function of the BFKL kernel, 
and 70 = 0.6275... solves x'ilo) = x(lo)/lo [7]- The front exhibits a universal 
tail: 

T(p, Y) ~ e -MP-Ps(Y)) for p _ ps > 1 (g) 
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These asymptotics set in diffusively and spread over the range p—p s (Y) within 
the time interval 

aAY ~ . (9) 

This process induces corrections to the asymptotic F-dependence of the sat- 
uration scale of the form 

dp,(Y) xho) 3 1 

— = " WoY ■ (10) 

Note that the velocity and the shape of the front for p ^> p s are completely 
determined by the linearized (BFKL) equation, and do not depend on the 
exact form of the nonlinearities: this is a very important consequence of the 
nature of the propagation of the front, which is pulled along by its tail. It 
implies that for a number of physical quantities, such as p s , we do not need 
to know the precise nonlinear mechanism that enforces unitarity of T. This 
means in particular, that the dipole picture is good enough for our purpose, 
although it is incomplete. 

Equation (10) shows that the velocity of the front increases with Y up to its 
asymptotic value, which corresponds to the velocity of a front which has the 
shape (8) all the way down to p — > oo. However, coming back to the physical 
situation where, due to the discreteness of the dipole number, the real front 
is a histogram, we see that T can assume the shape (8) only down to T ~ a 2 s . 
Starting from the initial condition at Y = and evolving it up to rapidity 
Y, the amplitude first grows until it reaches the unitarity limit T = 1 around 
r ~ r . Then the traveling wave front forms, and spreads from the point p s (Y) 
where T ~ 1 down to the point p at which T ~ a 2 . From Eq. (9) and from 
the shape of the asymptotic front Eq. (8), the latter process occurs within the 
rapidity interval 

r ln(l/^)/ 7 ol 2 



aAY 



(11) 



2x"(7o) 

(c is a number of order one) during which the velocity of the front keeps 
increasing according to Eq. (10). But once the point where T^a 2 is reached, 
the front cannot extend to even larger values of p (corresponding to lower 
values of T), at variance with the pure mean field case implemented by the 
BK equation. Accordingly, the front velocity cannot increase anymore. Thus 
the asymptotic F-dependence of the saturation scale is 

dp s (Y) _ x(7o) o _ 7oX"(7o) no x 

— ; = a oca — 7T- — — — . (12) 

dY 7o ln 2 (l/a s 2 ) 1 ^ 

The calculation of c requires a proper account of the exact shape of the front, 
and yields c = tt 2 /6 [14,22]. 

The result (12) is identical to the one obtained in the mean field approach of 
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Ref. [22], where the effects of discreteness have been simulated by introducing 
a second barrier on the phase-space for BFKL evolution (in addition to the 
first barrier meant to enforce saturation [10]), whose role looks indeed very 
similar to that of our above cutoff on the value of T. What was different, 
however, was the understanding of the physical role played by this barrier: 
in Ref. [22], the second barrier was merely intended to impose unitarity on 
those evolution paths which involve a single intermediate point. But it has 
not been realized there that, for the evolutions leading to a final amplitude 
Tf ~ 1 (as relevant for a study of saturation), the constraint implied by the 
second barrier is actually sufficient to guarantee unitarity for arbitrary paths 
(i.e. for paths with an arbitrary number of intermediate points). That is, the 
generality of Eq. (12) as being the correct result in the limit where Y — > oo 
and a s — > has not been recognized, nor argued, in Ref. [22]. 



4 The physical amplitude and its scaling 

So far, we have followed the evolution of the amplitude T(r) corresponding to 
one given partonic realization \u). Each such realization undergoes a stochastic 
evolution given by Eq. (6). At each rapidity, T(r) has the universal shape (8), 
up to fluctuations concentrated in its tail T ~ a 2 s . The position p s (Y) of 
the front exhibits the ^-dependence 5 given by Eq. (12). However, since the 
actual evolution is stochastic, the saturation scale p s undergoes a random 
walk about its mean (p s ), and p s gets a variance a. Physically, the origin 
of this phenomenon can be traced to the statistical fluctuations in the tail 
T ~ a 2 of the amplitude: in the course of the evolution, an unusually large 
number of dipoles may be created, which, after further rapidity evolution, 
would "pull" the whole front ahead of its typical evolution, resulting in a 
saturation scale p s for this particular realization larger than the mean (p s ). 
This mechanism leads to a spread of the saturation scales of different partonic 
realizations, while the shape of T in its bulk remains identical. In particular, 
the velocity of the average front d(p s )/dY is still given by Eq. (12) (at least, 
for sufficiently large Y; see below), because this is the (asymptotic) velocity of 
all the individual fronts making up the statistical ensemble. This is illustrated 
in Fig. 2, which in fact applies to the discrete statistical model in Ref. [14] 
(but a similar situation is expected in QCD): the profile of the amplitude 
for partonic realizations obtained from different stochastic evolutions over the 
same rapidity interval are shown. The dispersion of p s is manifest, as well as 
the universality of the shape of T. 

5 It is interesting to note in this context that, for a given front realization, the 
asymptotic velocity, Eq. (12), is reached exponentially fast in Y [17], at variance 
with the mean field problem, where the corresponding approach is only power-like, 
cf. Eq. (10). 
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The physical amplitude A(r, Y) is obtained by averaging T(r) over all Fock 
states at rapidity Y. In order to perform this average, the expression for the 
variance a of the saturation scale is needed. 



The latter was recently studied in the statistical physics context for various 
reaction-diffusion like models involving a finite number N of particles. The 
variance of the position of the front was seen to scale like 6 

- 2 = (Pi) - (Ps) 2 ~ (13) 
In (l/a 2 s ) 

from numerical simulations [14]. Although to our knowledge there is still no 
general analytical proof of this result and the status of (13) is still that of 
a conjecture, such behavior has been checked in independent numerical work 
(see e.g. Ref. [26]) and is also likely to be very general [17]. We also see on 
Eq. (13) that the fluctuations of p S) that are of order a oc V aY, are indeed 
subleading in Y with respect to the effects of discreteness discussed in Sec. 3 
(of order aY). This is a consistency check of the mean field calculation used 
to obtain (p s ), see Eq. (12) . 

Knowing the shape (8) of T, the value of the saturation scale p s and the 
amplitude of its fluctuations a, we are in a position to evaluate the physical 
scattering amplitude A(p, Y). Up to higher moments of the distribution of p s , 
A(p, Y) is obtained from the amplitudes T(p)\ Y for each particular realization 
of the Fock state of the target at rapidity Y (note that T(p)\y is implicitly 
a function of p s , as manifest e.g. in Eq. (8)), after averaging over the corre- 
sponding saturation momenta with a Gaussian weight of variance a (see also 
Fig. 2): 

A( ^=^J^ T ^^{ Jj1 ^) ■ (14) 
We deduce the following scaling form for the physical amplitude: 

Mp,Y)=A\ ;- ( "f» I , (15) 
Vay/ln 3 (l/al). 

up to 0(1/ VY) corrections. It is obvious from that formula that, at sufficiently 
high energies, geometric scaling does not hold for the physical amplitude. 
This feature is a direct consequence of the statistical nature of the parton 
model: the statistical fluctuations of the number of dipoles translate into a 



6 Formula (13) is borrowed from the third reference of [14], Eq. (4), with the re- 
placement cr 2 <-> Dn x t. The time t has to be replaced by aY and the number of 



particles N is the maximum number of dipoles of a given size 1/a 
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Fig. 2. The scattering amplitude T for different partonic realizations at a given 
rapidity against p = ln(rg/r 2 ). The thick line is the average over all realizations, 
i.e. the physical amplitude A, see Eq. (14). 

random wandering of the saturation scale, and after averaging over partonic 
realizations, the scaling form (15) results. 

The violation of geometric scaling manifest in Eq. (15) was already noted in 
Ref. [22], however, the square root in the denominator of the scaling variable 
was missing because of a lack of fluctuations in the tail of the distribution: 
the approach used there was relying on mean field throughout, missing the 
stochastic nature of the evolution. 



5 Conclusion 

We have shown that high energy QCD is similar to a reaction-diffusion prob- 
lem, well studied by statistical physicists. We have been able to obtain the In- 
dependence of the saturation scale, see Eq. (12), confirming results obtained 
recently by different methods [22]. We have also derived the scaling form of 
the asymptotic dipole-dipole scattering amplitude (15), which is related to 
the dispersion of saturation scales between different "events" (corresponding 
to different partonic realizations). That scaling is clearly not geometric. 

A recurrent theme in this Letter has been universality that guarantees that 
the lowest order results do not depend on the details of the model. The exact 
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way how saturation comes about was not an issue, as well the details of the 
elementary dipole interaction do not enter the leading order results (in Y ^> 1 
and C 1) that we have obtained here. Further terms in these expansions 
will be model dependent, and thus much more difficult to get. 

One of the points that remain to be studied is how fast the computed asymp- 
totics set in. The BK equation may still be a good approximation for a large 
target (like a nucleus) and in the first stages of the evolution, when the trav- 
eling wave front has not diffused down to T ~ a 2 s . More precise numerical 
studies [18,21] may help to clarify this point. 

Finally, as mentioned in the Introduction, the dependence upon the impact 
parameter plays an important role for the overall physical picture of unita- 
rization. This will be discussed at length somewhere else [23]. 
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